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^ Abstract 

▼— H Interacting massive fields with m > -^ va. d + 1 dimensional de Sitter space are fundamen- 

1^ tally unstable. Scalar fields in this mass range can decay to themselves. This process (which is 

*^ kinematically forbidden in Minkowski space) can lead to an important change to the propagator 

%^ and the physics of these fields. We compute this decay rate by doing a 1-loop computation for a 

massive scalar field with a cubic interaction. We resum the 1-loop result by consistently solving 

the Schwinger-Dyson equations. We also perform an explicit resummation of all chain graphs in 

the case of the retarded propagator. The decay rate is exponentially suppressed for large m/H 

and the flat space answer (vanishing decay rate) is reproduced in that limit. 

PACS numbers: 



I. INTRODUCTION 

Much work has been done in studying the quantum behavior of low mass {m ^ H) 
scalar fields in de Sitter space (dS). In the infiationary paradigm, quantum perturbations 
of these scalar fields give rise to the large scale structure that we see in the Universe today. 
Comparatively less attention has been paid to massive fields with m ^ H. Clearly in 
the limit of very large mass or vanishing Hubble rate, one should recover the fiat space 
answer. Recently, however, Krotov and Polyakov pQ analyzed the correlation functions of 
these massive fields and found a surprising result. They found that in the Poincare patch, 
the Keldysh propagator F at 1-loop receives a 1-loop logarithmic correction 

F(fc, r) -^ F^{k, r) (1 + cT In kr) , (1) 

where F^{k, r) is the tree level propagator, k is the spatial momenta and r is the conformal 
time. The loop correction naively diverges at late times r — ?► 0^, which naively suggests a 
breakdown of perturbation theory. 

This 1-loop diagram had been analyzed before by Marolf and Morrison [2] working on 
the sphere and analytically continuing to obtain the result in de Sitter space. By resuming 
all the IPI diagrams they argued that this logarithmic correction actually corresponds to 
an imaginary shift to the mass 

fi ^ fi + ia (2) 

where /^^ = ^ — ^ > where H is the Hubble constant and d is the number of spatial 
dimensions. In real space, the main effect of the resummation is to make the propagator 
decay faster at large distance 

GiZ)^^^{AZ"^ + BZ-^n, (3) 

where Z is the geodesic distance. A, B are specific coefficients and a will turn out to be real 
and negative. 

In this paper we perform a resummation of the diagrams directly in Lorentzian signature. 
We find that it is possible to resum all the bubble diagrams to show that the 1-loop answer 
indeed leads to a (time independent) imaginary shift of the mass. The resummation is 
performed by consistently solving the Schwinger-Dyson equations for the system using an 
ansatz for the full resummed propagators. In the case of the retarded propagator we directly 
carry out the resummation of all bubble diagrams and show that the result is in agreement 
with the solution of the Schwinger-Dyson equations. We also argue that the same imaginary 
contribution can be found for quartic interactions. For low mass fields, very similar logs 
appear and these have already been resummed using a variety of techniques (for recent 
work see [3H9]). Our computations are all done in the Poincare patch of de Sitter and our 
resummation agrees with the IPI resummation on the sphere [2] confirming that the two 
computations are equivalent [TO] . 



We note that this phenomenon may be related to another fact: interacting massive fields 
in the principal series (m > ^) appear to be fundamentally unstable. Since conservation 
of energy does not hold in dS, it is possible for a particle to duplicate itself or even decay 
into heavier ones [TT]. Interestingly this is only true when m > ^. 

To our knowledge the computation of this decay rate was first done by Bros, Esptein and 
Moschella in a series of papers [111 [13] (see also (THUS]). In these papers, they computed the 
tree level decay rate for a scalar field of mass m with cubic interaction A0^. In flat space, the 
decay rate is exactly zero simply because of kinematics. In de Sitter space, however, they 
found that there is a non-zero but exponentially small decay rate (for the principal series). 
This is a reflection of the energy non-conservation which opens up new decay channels. In 
our case, we have found an imaginary shift in the mass, and it is plausible to relate this to 
the decay rate found by [121 [IS] • 

We compute explicitly the decay rate (aka imaginary part to self-energy) and we show 
that it goes to zero roughly as A^ (^) q-'^"^/^ in the large mass, low H limit where the 
exponent n depends on the number of dimensions. The decay rate goes to zero very fast in 
the large m/H limit and the flat space answer is recovered. The decay of the Higgs (say 
with m ~ lOOGeV) in our Universe today (with Hq = 10~^^ev) is obviously very small. For 
m r^ H the decay rate reaches a flnite value of order F ~ A^ . 

This paper is organised as follows. We set up the problem in section § [IT] and compute the 



1-loop contribution in section § |III| for all propagators. We then proceed to resum using two 
different methods in section § IV Finally, the decay rate is explicitly evaluated in section 
§ |V| In an appendix at the end, we show how a fleld with quartic interaction will have similar 
loop corrections. 

II. INTERACTING SCALAR FIELD THEORY IN DE SITTER SPACE 

We are interested in calculating in-in correlation functions. 

(Q(t)) = {in\ U^QU \in) (4) 

where the evolution operators are 

[/(t,to) = T exp (-i I dt'Hi{t')\ , U\t,to) = T exp U f dt'Hi{t')\ . (5) 

We take to = ti„(l =F ie) where the sign must be chosen accordingly to select free vacua in 
the far past and Hj is the Hamiltonian in the interaction picture. One can compute the 
1-loop contribution by directly expanding the evolution operator to the desired order. For 
the purpose of resummation, we flnd it useful to use Feynman rules that can be obtained by 
doubling the number of flelds in the Keldysh-Schwinger formalism. For a scalar fleld with 



non-derivative interactions of the form 

H,{t)= d'x^gY,-^. (6) 

n 

the Feynman rules are well known. In the context of cosmology where the metric depends 
on time, it is useful to have either a full real space representation or a mixed representation 
where we use momentum space for the spatial direction but keep the time dependence 
explicit. (For more details on how to find the Feynman rules in a cosmological context see 
for example [12] or the appendix of [H]). We will use the Keldysh basis for the propagators 

G^'ix.y) = 9{xo - yo) (G-+(x,y) - G+~{x,y)) , 

G^{x,y) = 9iyo-xo){G^-ix,y)-G-^{x,y)) , (7) 

where F stands for the Keldysh propagator and G^'"^ are the retarded and the 
advanced propagators respectively. The Wightman functions are G^^{x,y) = 
i{(f){x) (f){y)) ,G^^{x,y) = i{(l){y) (j){x)). In our computation, we will use exclusively the 
mixed representation where we Fourier transform only the spatial part of the propagator. 
We show the resulting Feynman rules in Fig. ([I]). 

The propagators can be obtained from the mode functions 

0(^, t) = I ^ {a,e^'^-%{t) + a}/^-^\j; (t)) . (8) 

We will work in the Poincare patch of dS with metric 

ds = (9) 

where we have set the Hubble constant equal to 1 and the conformal time r runs from — oo 
to 0. The number of space-time dimensions is -D = (i + 1 . In the Bunch-Davies (BD) 
vacuum, the modes are jT] 

„-„i7r/4 nzi _-.\dl2 i 

f/.(r) = %^ ^ e-^/^<^(a;) ^ i=(-r)^/^Mx,/.) , (10) 



where x = —kr, /x = ^Jrr^^^^{(ipX^ , we have dropped the irrelevant phase and 

M^,/^) = y|e--^/2i7«(x). (11) 

We work in the principal series where [i is real and the order of the Hankel function is 
imaginary. The function h is symmetric in the mass h[x^ /i) = /i(x, — /i) and has the following 
asymptotic behavior for large and small values of x 



ie*(a:-f) -X » 1 

h[x,y.)^{ V- ' (12) 

Aj.x'^ + y4_a;"*^ ; a; < 1 
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FIG. 1: Propagators and some vertices (A'^ = 3 and A'^ = 4) for a scalar with a potential of 
the form XIat ^vT'^^- ^ote that Gyi(/c,ri,r2) = Gji{k,T2,Ti). We use the shorthand notation 
^12 = ^(''"1 — T2). When a two-point function is attached to a vertex, the corresponding time must 
be integrated over. Internal spatial momenta must also be integrated over f .^ ?^ . As usual, one 
must also divide by the symmetry factor for the diagram. At A^ = 5, there are 3 vertices, etc. The 
scale factor in dS with iJ = 1 is just a{T) = —1/t. 



where A± are specific functions satisfying the following useful identities 

|A+|2-|A_|2 = i, |A+|2 + |A_|2 = icoth(7r/i). (13) 

In terms of those modes, the propagators are 

F{k,n,T2) = {nT2Y'^f{kr,,kr2) = ^{nr^Y'^ ^ihi-kn,fx)h*i-kT2,fx)) , (14) 

G'«(A;,ri,r2) = 9in - r2)r^/\^/'g''ikn,kT2) = -{nr^Y^'^ein - r^) '^ ih{-kn, fi)h*{-kT2, fi)) . 

We will often factor out {tiT2Y^'^ as well as any time ordering function and denote the rest 
of the propagator by non-capitalized letters g^ and /. We can expand the propagator for 



small value of the arguments — A;ri^2 ^ 1 using Eq. (12) to get 



/ ^ ^ l2A+A*_{k'nT2r' + 2A_A\{k\r2)-'^' + - coth(7r/i) ( ("^^ '^ + (7^ 

We see that for small kTi^2 both propagators have k independent oscillating terms. In the 
loop computation these oscillating terms can interfere destructively. When this happens one 
gets a secularly divergent log in time Infer. We will compute this log contribution to the 
1-loop diagrams below and we will then proceed to resum the 1-loop bubble diagrams. 

III. THE CUBIC 1-LOOP CORRECTION 

In this section we compute the 1-loop correction from a cubic interaction to the retarded 
and the Keldysh propagator. 

A. The Retarded Propagator 

Let us start with the 1-loop correction to the retarded propagator from a cubic coupling 
-^. There is only one diagram which is shown in Fig. (2). The 1-loop correction to the 




T2 

FIG. 2: The 1-loop (f)^ correction to the retarded propagator; this is —iG^^{q,Ti,T2). 

retarded propagator is 

G"(,.r.r.) = A^/^/l*3/^r.^G«°(,.r,.r.) (16) 

xG™(*:,rs,T4)F"(|i-?l,T3,T4)G™(<(,T4,T2), 

where G™ is the tree level propagator while G^^ is the 1-loop part. A similar notation 
is used for F. There is no log coming from k = since the field is massive. Instead, a 
logarithm arises from the region where k ^ q. To see this, we first factor out the overall 



factor {tiT2Y^'^ from (16) to obtain 



/°° d'^k l'^'^^ f^^ a 

— — ^ / rfX3 / rfX4(x3X4)'/'-'(7™(gri, |X3)(?™(X3, X4) 

x/°(a;3,X4)^™(|x4,gr2). 

where we have defined x^ = kr^ and X4 = A;r4. We expand the two external propagators for 
~Q"'"i,2 ^ 1 and —f 3^3,4 ^ 1 as 

In the product of the two (7™ above, there are two terms which interfere destructively 



.2/i/ Jq (27r) k J^^^ Ji^^^ 



g-' ^ XU \ / — — — / C?X3 / dx4ix3XiY/^ ^g^%X3,Xi)f{x3,Xi) 



In the limit that —kri — )■ and —kT2 — ?■ cxd, the time integral is independent of k 

fcn px-i pO i-xj, ^ pO pO 

dxs / dxi ~ / dxs / dxi = ^ dx^ / (ix4 . (19) 

kT2 J kT2 J — 00 J —00 ^ J —00 J —00 

Note that this is only valid for —kri -^ 1 and —kr2 ^ 1 which imposes new limits on the 
k integral. Since (7™ is odd under the exchange of two arguments while /° is even, the two 
terms are the same up to a minus sign. Hence 



9"'-^ ^ -P) 1/ ^. (20) 



ifi /_ \ JM\ />- l/ri ^^ 

2/Lt \^V^2y \rij J J-i/r2 k 

with 

iX^jaf{^i)Sd-i 



a = 



pO 



.Pf(^) = 1 / dx, I dx3X~'^xTix3X,Y/^-'g''%X3,x,)f{xs,x,), (21) 

where Sd-i is the volume of the d — 1 dimensional sphere. As promised this part of the 
diagrams is UV convergent and the momenta integral only give rises to a logarithm of ratios 
between the two times. 



Altogether, the retarded propagator at 1-loop in the hmit ti > T2 and for — gn 2 ^ 1 
goes hke 

G^(°+^)(g,ri,r2) ^G™(g,ri,r2) [l-aln^l . (22) 

L ''"2 J 

Some comments are in order regarding the (multiple) approximations done in obtaining this 
simple answer. We evaluated this loop contribution assuming ti ^ T2 and we took the 
external momenta to be small with q <^ — 1/ri. We neglected the k < q part of the integral 
as well as approximating the limits in the J{fi) integral. Both approximations will lead to 
extra k dependence which then lead to power law contribution of the type 0{ti/t2) (which 



are negligible in the late ri — )■ limit). Using Eq. (14) 
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16 ./o 



J^''^(;u) = — / dx4dx3X^^^x^l'{xsX4)'^''^ ^{h{x3)h*{x4) — h{xs)*h{xi)) 
x{h{x3)h*{x4) + h{x3)*h{x4)) 



with 



= -i^{\9if^)\'-\9i-f^)n, 

= ^l^(-/^)P(l-e-^-^), (23) 

gifx)= dxx'^/^-^+'''h{x,fx)\ (24) 

Jo 

This is the same formula as in [1]. The last line was obtained using \g{fi)\'^ = e~'^^^\g{—fi)\'^ 

which follows directly from the properties of Hankel functions. The 1-loop coefficient is then 

given by 

^-^M-^>^'^^-^-'n, (25) 

which is real and negative. In section § IVl we will analyze this integral further. 



B. The Keldysh Propagator 

There are three 1-loop graphs contributing to the Keldysh propagator as shown in Fig. ([s]) 
(there are five diagrams including time reversals). These diagrams may be computed using 
the same approximations and methods as in the last section. 



For all diagrams, one obtains equations similar to Eq. (18) with the following differences. 
First, we have either two g^, two / or one g^ and one / in the time integral. We denote 
the new integrals as J^^ and J^^^ respectively. Second, the limits on the time integrals are 
different. For diagrams (1), the limits of integration are 

dxs / dx4 . (26) 





I) 



2) 



3) 



FIG. 3: Corrections to the F propagator. The time reversal (reflection) of all diagrams needs to 
be included as well. 



and in the time reversal, we replace ti by T2. For diagram (2), the limits are now 

dx^ / dx4, . (27) 

—00 J —00 

and this is unchanged for the time reversed diagram (we took ri > T2). Finally for diagram 
(3), the limits of integration are 

fcTi rkT2 

dxs / dx4 . (28) 

■00 J —00 

and there is no time reversal. Demanding that these integral be fc-independent require both 
—kTi2 ^ 1. This gives an upper limit in the momenta k but the lower limit remains at q 
(to be contrasted to what happen for the retarded propagator). 



The Keldysh propagator in the small — gr ^ 1 is given in Eq. (15); here we further 
neglect the qr pieces compared to the other terms 

/Vn,r,)=^^f(^V%(^V*''V (29) 



4/i \ V^2/ \T2 

Approximating /° to this term, one finds that the 1-loop correction from diagram (2) and 
(3) is 

/' = -f{q. n, r2)(a2 + as) In -qr^ , (30) 

where again ti > T2 and 

^^ = 4;.coth(7r'/^)(2vr)^^''^^^' ^^ = -;.coth(/^K2vr)'^^''^^^ ' ^^^^ 

We would like to emphasize the relative sign between a2 and a^ due to the presence of 4 
retarded propagators (each with their —i factors) as opposed to 2 and also note the factor of 



1/4 from the vertex. Also we have included the time reversed diagram for 02 which amounts 
to an extra factor of 2 and there is a symmetry factor of 1/2 for both diagrams. Now both 
(72 and (73 are real and they lead to an imaginary shift to the mass. It is worth pointing out 
that the sum of the two terms does not contain a divergence of the type J h{x)h*{x) due to 
the identity 

Uff-J'' = \{\gifi)\' + \gi-li)n. (32) 

Finally, diagram 1 in Fig. ([s]) also has a very similar structure. We expand the Keldysh 
propagator in the last leg as before and the final result is 

/(°+i)(g,ri,r2) = f{q,n,r2) [1 - M\n{-qn) + H-qr^))] , (33) 

with 

zAV^/(/x)g,,i 
^^ = (2vr)'^2;. ' (^^) 

Adding all the contributions together and writing ln(— gri) = In — + 21n(— gr2), we get 



/(°+i)(g,ri,r2) = r(g,ri,r2) 



1 - o"! In (2(Ti + (72 + (73) ln(-gr2)) 

T2 



(35) 



Using the properties of the integrals J^^{fi), J^^{ij) and J^^^jj) in Eq. (23) and Eq. (32), 
one can show that 

2ai + (72 + ^3 = . (36) 

Therefore the full answer is just 

f^'^'\q, ri, T2) = f{q, n, r^) [l - o. In ^"j (37) 

and (7i is the imaginary shift to the mass. It is exactly equal to the shift we found for the 



retarded propagator in Eq. (21) 



IV. RESUMMATION OF LARGE LOGARITHMS 

In the last section we computed the 1-loop contribution to the retarded and to the Keldysh 
propagators in A0^ theory. We found that the result contains a potentially large logarithm 
that is either a ratio of two times or a ratio of the latest time to the external momenta. 
In order to obtain this result we had to assume that ri ^ T2 or that gri^2 ^ 1 and the 
logarithmic contributions are then dominant and large. These large logs cause trouble as 
they naively signal a breakdown of perturbation theory at late time. To address this problem, 
we will now resum these 1-loop logarithms by using two different methods, both of which 
give the same result which is finite and well behaved at late times ri — )■ 0~. 
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The resuinination technique demonstrates that the 1-loop contribution resums to a time 
independent (in this case imaginary) shift to the mass. Indeed, if one performs this com- 
putation by Wick rotating to the sphere as in |2j, there is no logarithmic time dependence 
and the 1-loop simply gives a finite contribution to the self-energy which shifts the mass. 
So the results of our resummation shows clearly that the same answer is obtained using 
the in-in formalism in the Poincare patch versus using the sphere. On the other hand, our 
resummation methods can easily be extended to non-equilibrium situations that may not 
be captured by a computation on the sphere; for example global de Sitter space [T] or an 
unstable vacuum (as opposed to the Bunch-Davies vacuum). 

In the case of the retarded propagator we are in a fortunate situation because there is only 
one diagram at 1-loop and therefore it is a simple matter to resum. Large infrared logarithms 
have been argued to be well captured by the classical approximation (see for example [16]) 
or by semi-classical methods such as the stochastic approach [18] (see also [T91 - I2T] ). Here 
are two different ways to resum the 1-loop correction to the retarded propagator where the 
errors are easily quantifiable and the next order correction is in principle calculable. 

1. In the first method we will set-up a Schwinger-Dyson like non-perturbative equation 
and solve it by making an ansatz for the 1-loop exact retarded propagator. 

2. In the second method we explicitly compute the n— loop bubble diagrams and we 
resum them all (this is for the retarded propagator). 



It would also be interesting to use dynamical renormalization group methods (see 
and [HI [23 ES] for some applications in cosmology) to perform the resummation but we leave 
this to further work. 

A. Schwinger-Dyson Equation for the Retarded Propagator 

Considering only the bubble diagram, it is simple to set-up a non-perturbative equation 
for the fully resummed propagator a la Schwinger-Dyson (see Fig. |4]). The resulting equation 



is of the form Eq. (16) with the last propagator replaced by the full resummed one. We 



make the ansatz that the mode functions, which were originally written as Eq. (10), should 
become 

U,iT)^{-TY/'h^ix,fi), (38) 

where h^{kT) are some new functions. Correspondingly, the propagator is now of the form 

G^{q,n,n) = e{n - T,){nT,Y/'g^{qn,qn) . (39) 

We can write the Schwinger-Dyson equation as 
g^{qTi,qT2) = g^\qTi,qT2) 

^ y dT,j dnir.nY^'-'g^^qn, qT,)g^'ikT„ kT,)f{kT,, kn)g^iqu, qr^) , 

11 
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FIG. 4: The Schwinger-Dyson equation for the retarded propagator. The hne with the square is 
meant to represent the fully resummed propagator. 



where as before we have taken the Umit k ^ q. We perform the Fourier expansion of 
g^{qT4, qT2) in terms of 5x4, and in particular, we consider a Fourier component which is of 
the form {qTiY°'ha{qT2) ■ The modes a = ±fi received a large enhancement from the integral 



proportional to lnri/r2 as we found in Eq. (22). Now ex hypothesis the resummed propagator 



should have a good perturbative expansion for all time, and hence we conclude that there 
are no Fourier components of g^{qT4,qT2) proportional to (5x4)=''*^. Since at leading order, 
such terms are indeed present, what must therefore happen is that the exponent of these 
Fourier components is corrected by a shift which is proportional to the coupling. 

The mode functions which receive the largest corrections are thus the ones proportional 
to (f){k, t) ~ (A;r)^*^; we assume that these mode functions are corrected to be of the form 
0(/i;, r) ~ (kr)'^^^'^ for some b. These will produce large corrections to the propagator which 
depend on the integrals 

I±= I TTZ-J dxJ dx^{x,x^Y'^-'xf^xf^g^\x,,x,)f\x,,x,)xt 

) "' J —00 J —00 



where J^-^{±jj,b) is the integral defined in Eq. (21) with the extra xj in the integrand. 
We have employed the same approximations as in section § III and ti ^ T2. We split the 
resummed propagator into the resonantly enhanced and the resonantly unenhanced modes 
by 

A,U.,r.)-±(^{^J-[^y^ (^)".-,«fen.,.,) (41) 

and 7"^ contains no terms with Fourier component (5x4)^*'^. The terms in 7^^ does not lead 
to large integrals on the right hand side, and can be dropped at leading order on the RHS. 

12 



Since every other term in the equation only contains the Fourier components r^ *^, it follows 
that 7^ is zero at leading order. With this ansatz for the full retarded propagator, the 
Schwinger-Dyson equation becomes 



-.R ^ „-Ro _|_ m ^(") 



(42) 



which of course exactly reproduces the 1-loop answer Eq. (22) when 5 = 0. The factor a(b) 
is as given in Eq. (21) with again the extra factor of X4 in the integrand 



a{b) 



,^ w.^ / dX4 / dX3X^'^xf{x3X4Y^'^''^g^^{x3,X4)f{x3,X4)x^^. 

(zvrj 4/i j_oo J -00 



Now we see that we can consistently solve the SD equation by setting 

b = a{b) 



(43) 



(44) 



Which in general is some transcendental equation to be solved for b. We can approximate 



a{b) ~ a (given by Eq. (21)). This amount to resumming only the bubble diagrams and 



nothing else. The upshot of this computation is that the ansatz Eq. (41) consistently solves 



the SD equation for b = cr(6). This modification of the propagator can be seen as a modifi- 
cation of the modes. As a result, the same ansatz can be used for getting the resummed F 



propagator. We will rederive this result in subsection § |IVC[ by doing direct resummation 
of 1-loop chain diagrams for the retarded propagator. 



B. The Keldysh Propagator 

We now turn to the F propagator. Unlike the retarded propagator, we have more bubbles 
to account for in the Schwinger-Dyson equation, see Fig. ([s]). The structure is completely 
similar as in the G^ case. For example, including only the first diagram, the SD equation 
reads 



/(?n,gT-2) = f{qTi,qT2) 

j^,j dnj du{Tsnr/'-'g''\qn, qT,)g''\kT,, ku)fikn, kn)f{qu, qr^) + 

We can consistently solve the SD equation by starting with the same ansatz as we had for 
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R 



fiqn, qT2) 



coth(7r/i) ( ( T4 



4/i 



T"2 



ifi 



+ 



-?./i' 



-b 



(45) 
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FIG. 5: The Schwinger-Dyson equation for the Keldysh propagator 



for some undetermined b. There is five 1-loop computations to perform and we denote them 
1 to 5 going from top to bottom in Fig. (|5| 



/i = a^ib) 



f 



-qr2Y 



(46) 



J\ 



/3 = o^3(&)y(l-(-gr2)^): 

u = -^ — r(l- (-9T-2) ) 



h 



2 b 

<y2{b) f 



(1 - i-qn)') . 



Each crj(&) are exactly as given in section (IIIB) but with an extra factor of x^ in the 
integrand. The SD equation then becomes 



f = r + f 



, .O^lW / /T-2 



n 



1 



f 



(2(7i(6) + (72(6) + (73(6)) (l-l-grs)") 



(47) 



where the first term is the tree level contribution and we have rearranged the contribution 
from the first and second diagram. As before one can show that for any b 



2(Ti(6)+(T2(6)+(T3(6) = 



(48) 
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The SD equation can consistently be solved by setting 

b = ai{b). (49) 

Neglecting all higher order terms and keeping only bubble diagrams we have that ci = b. 
Note that the shift to F is the same as the shift to G^, that is ai = a. 

C. Direct Resummation 

A direct computation of the 2-loop chain diagrams for the retarded propagator is simple 
to do, one simply replaces the last propagator by the 1-loop result 

oo Jdh. rkri pxa 



-,R2 



r^" = ^'/ TTT^TTTw / dxs dx^{x^x^Y'^-^g''\qT^,}x^)g''\x^,x^) 



x/°(x3,X4)(7^i(|a:4,gr2). (50) 

The 1-loop propagator contains a lnA;r2 which when integrated over k will give a (In— )^. 
In the large time difference limit |Ti/r2| ^ 1, the log term will dominate over other 
corrections (it is not enough to simply demand —qTi^2 ^ !)• Finally we approximate 
In j^ ~ — ln(— A;r2). The lnx4 is dropped because it contributes a subleading lnri/r2 in 
the final answer. Hence 

-l/ri ^dj^ 

(51) 




2 
A clear pattern emerges and the n— loop chain (bubble) diagrams is just 

g^- = g^'a-^^, (52) 

n! 



and the sum exponentiates to give 



gR ^ ^i?Og-a In n/ra 



2^ 



ILL— (J / \ —til — a 



T2 \T2 



(53) 



The direct resummation of the Keldysh propagator does not work as well partly because of 
the gr terms in the propagator. This is not really that surprising as it often happens that 
one can consistently solve the SD equations (as we did for the Keldysh propagator) but it 
still proves impossible to perturbatively resum diagrams. 
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V. THE DECAY RATE 



The resummed propagator (for e.g. the retarded propagator Eq. (53)) has exactly the 



same form as the bare one. Using mode functions representation from Eq. ( 14 ) 

^^(fc,ri,r2) = '^{h{-kr^,^^)h\-kr,,^^)) = "^Q (<L(^)4^-.).(^)) (54) 

we see that the resummation of bubble diagrams has lead to an imaginary shift to the mass. 

fi ^ fi + ia , (55) 

with 

'-i^M-^n^-^-n^-'i^K(,). (56) 

Note the complex conjugation of the order parameter in the second Hankel function; one 
should not simply shift the mass in the propagator since there is complex conjugation going 
behind the scene. This formula is valid for any /i real and greater than 0. 

li ^ ^ m ^ 0, then the propagators approach their flat space equivalent and the imagi- 
nary shift to the mass leads to a Breit-Wigner resonance of the form 

^ (57) 



p2 _ j^2 _|_ ^.jj^Y 

Therefore, in analogy with flat space results, we define the decay rate to be 

r = -2cT (58) 

which is positive (as it should be) since a is negative. It is simple to perform the integral g{n) 



given by Eq. ( 24 ) numerically. The integrand oscillates wildly for large x and we regulate it 
with an exponential cutoff. In Fig. ([6]), we plot the function K(fx) as a function of fi. We 
find that the decay rate is Boltzman suppressed as found in [13]. 



ir(/i) > e-""" . (59) 

We find that the integral converge in any dimensions. For fi = oi m = ^, the decay rate 
is finite and of order 

r«0.05^. (60) 

The angular factor of Sa-i/i^nY is as usual for loop computations. 

VI. CONCLUSION 

We have shown that there is indeed an imaginary shift to the mass of a massive field with 
a (j)^ interaction, as argued previously by Marolf and Morrison ^. In the in-in Poincare 
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FIG. 6: The integral K{fj,) as a function of // for d = 3. 



calculation, the loop contribution receives an imaginary contribution that diverges at late 
time. We have shown that one can consistently solve the Schwinger-Dyson equation by shift- 
ing the mass by an imaginary component. This effectively resums all the bubble diagrams 
directly in the Lorentzian signature spacetime and reproduces the IPI resummation from 
the sphere. 

Many other resummation techniques could be used and it is worthwhile to explore them. 
From the Schwinger-Dyson equations, one can set-up a Boltzman kinetic equation. Alter- 
natively, stochastic techniques [ISj or dynamical renormalization group techniques [HI [25] 
could also prove to be very useful. 

Our results give additional support to the statement that interacting massive fields with 
m > ^ in de Sitter space are fundamentally unstable to decay. The reader may worry 
about the interpretation of this imaginary shift to the mass as a decay rate given that in de 
Sitter, there is no well defined notion of particles and the S-matrix is an ill defined quantity 
(the Jost vacua of [2Z] - see also |2H]^ is the exception to this statement). In fact there is no 
contradiction as a decay rate is not an in/out quantity. One can never take a fundamentally 
unstable particle to t — t- — oo to prepare an initial state. An unstable particle is not really 
a particle, it is a resonance. If one take the large mass limit of the resummed propagator in 
dS one recovers the flat space answer 

1 



pz 



rn^ + imV 



(61) 



and clearly our imaginary shift to the mass then lead to the usual Breit-Wigner resonance. 
Our notion of decay rate differs from the one used in |25l EHl ED] where they are looking 
for time dependence in correlation functions (in particular, they are looking for decay of the 
vacua to some other state). The case of conformally coupled scalar with cubic interactions 
was recently examined in [31] . In that paper, they both found a time dependent anomalous 
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dimensions and a decay rate. For the quartic potential, they find no such decay rate. Again 
to contrast, we have examined the case of a massive field in the principal series which is 
perturbatively stable for cubic or quartic interactions but both cases show an imaginary 
contribution to the self-energy which we interpret as a decay rate. 

The case of global de Sitter is very interesting. Krotov and Polyakov also computed the 
1-loop diagram in global de Sitter. They argued that the euclidean theory does not correctly 
capture the physics in this case. As global de Sitter can roughly be viewed as two copies of 
Poincare, the result for the 1-loop is very similar 

F(fc,r) ^F(A;,r)° M +aln^ j (62) 

where tq is the time at which one turns on the interaction in the lower (contracting) Poincare 
patch. This result was obtained in the approximation that r ^ tq. This loop contribution 
does not seem to simply resum to a shift in the mass and may be a signal of more interesting 
dynamics. In fact, Krotov and Polyakov argue for a decay of vacuum energy along the line 
of [271 132] (see also [331 El] and [35ti37] for similar claims). We however note that the result 
of [T] is only valid for r ^ Tq which appears to imply that there is no large logarithm in 
this case. Recently, Marolf and Morrison [38] have examined the question of what happens 
at late time and they found that the full propagator approaches the Euclidean one. This 
is in agreement with the proof of IR stability on the sphere to all order by Hollands [3S] 
and Marolf and Morrison [30] (see [H] for a discussion of the massless case). It would be 
interesting to see if our methods could be apphed to this case to find a better result for the 
propagator. 

Finally, we gave an explicit formula for the decay rate. It is suppressed by the Boltzman 
factor at large mass F — )■ \^q-'^^/h but it reaches a finite value of order F ~ A^ when 
m = ^. This is in agreement with the tree level results of [12]. It is tempting to interpret 
this decay as coming from the "thermal" nature of de Sitter. In this interpretation, fields 
with m > d/2H have a Compton wavelength smaller than a Hubble patch and one can use 
the static coordinate to describe them. In this system of coordinates, the Bunch-Davies 
vacua looks thermal and one can naively steal energy from the heat bath to open new decay 
channels. This intuition is not borne out by analog computation in thermal field theory (see 
|42j ) . While there is in general an imaginary part to the self-energy at finite temperature 
and energy can be taken from heat bath, the process of a massive fields with m ^ T 
self-duplicating themselves does not occur. 

We reiterate that this imaginary part only shows up when the mass goes above the upper 
bound 

m^ < ^ . (63) 

Below this bound the self-energy is real and there is no decay rate (see for example [H]). This 
bound is of course very similar to the well-known Breitenlohner-Freedman [ISl HI] bound 
in Anti-de Sitter 'ni?R\^g < — ^. The de Sitter analog was proposed by Skenderis and 
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Townsend |l5] using arguments from fake supersymmetry. When one attempts to define a 
holographic dS/CFT correspondence (see for example [IHIIIT] and references therein), fields 
with mass above this bound are the "harder" case to interpret. The conformal weight of 
these fields living on the boundary are imaginary pointing to a non-unitary CFT. Here, we 
see that when interactions are turned on, the fields now decay even faster at large distances. 
It would be interesting to work out the implications of this decay rate for dS/CFT. 

Finally, beside these very theoretical questions, this decay rate could have practical im- 
plications for the phenomenology of inflationary theories. Fields with m ~ if could be 
present during inflation and they could have an effect on the spectrum of perturbations. For 
example, in the model of quasi-single field inflation |18] , oscillations of a massive field lead 
to features in the power spectrum. In [48.J, the authors restricted themselves to m ~ if but 



still less than |ii and it would be interesting to see what changes if the mass is above this 
bound. Another change to the physics may occur when integrating out massive fields to get 
the effective action at the inflationary scale. This is usually done as in flat space, or at the 
very least neglecting interactions (see [M] for a recent analysis). Integrating out massive 
fields with this new modified propagator may affect the quantum corrections in the effective 
action. We hope to return to these questions in future work. 
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Appendix A: Sunset diagrams in X^cj)'^ 

The imaginary contribution to the self-energy is not peculiar to the cubic interaction. 
Here we show that for a quartic interaction of the form A40^, there is an imaginary component 
as well. (The subscript is just to differentiate the quartic coupling from the cubic coupling 
used in the main text). The relevant diagram is the one depicted in Fig. ([T]). This is the 
only sunset diagram at 2-loop for the retarded propagator. It has a very similar form to the 
cubic case 

g^' ^ Xl rd%d% r ^ r '^-^{r,T,Yg''\qn,qr,)g^\k[T,Mu) 

Jq Jt2 ^3 Jt2 ^4 

x/0(A;^r3, hn)f\\k[ + k',\Ts, \k[ + k',\T,)g^\qn, qr^) , 

where we have taking fci 2 ^ g as before. Rescaling times in terms of fci 2 as we did in section 
is not particularly useful here because of the norm \k'i + /C2I. Instead, it is easier to 
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FIG. 7: 2-loop correction to the retarded propagator in Xcj)^. 
rescale the momenta 

We expand the two external retarded Green's function and we keep only the interference 
term which contain ratios of T^/Ti- As before the remaining time integral is antisymmetric 
under the exchange of time and so the end result is 

.- - Ai f^v ii^^r - f^n r ^ r ^ r^%<^% m 



where we use a shortand notation g^^{ki^T^/Ti) = g^^ ( ^ia/?' ^ia/t^ )• "^^^ integral over 
fci 2 will lead to a function of T3/r4 only; call this function $. The integral over time is 
logarithmically divergent 

"^r^*,,/..)^rfP^*W^.n-. (A3) 

2 ^3 Jra ^4 Jt2 ^ Jl Z T2 

Finally the integral over momenta is imaginary because g^ is imaginary while / is real. In 
the end, the sunset diagram goes like 

^^i~^«°A^^/(/z)ln^, (A4) 

/i T2 

where /(/x) is imaginary and equal to 

/Ti/r2 J /"OO 

-z'^ j d%d%g''\h, z)fik2, z)fi\h + k2\,z). (A5) 

The 2-loop sunset diagrams in A40^ has exactly the same structure has the cubic diagram 
of Fig. ([2]). We can solve the Schwinger-Dyson equation consistently using exactly the same 
ansatz as we did and the end result is an imaginary shift to the mass. This shows that the 
decay rate is not particular to the cubic case with its potential unbounded from below. 
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